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Abstract 



Higher order Delaunay triangulations are a generalization of the Delaunay triangulation which provides a 
class of well-shaped triangulations, over which extra criteria can be optimized. A triangulation is order-fc 
Delaunay if the circumcircle of each triangle of the triangulation contains at most k points. In this paper 
we study lower and upper bounds on the number of higher order Delaunay triangulations, as well as 
their expected number for randomly distributed points. We show that arbitrarily large point sets can 
have a single higher order Delaunay triangulation, even for large orders, whereas for first order Delaunay 
triangulations, the maximum number is 2™ . Next we show that uniformly distributed points have an 
expected number of at least 2 P1 "' 1+ °( 1 " first order Delaunay triangulations, where pi is an analytically 
defined constant (pi ~ 0.525785), and for fc > 1, the expected number of order-fc Delaunay triangulations 
(which are not order-i for any i < fc) is at least 2 Pfcn< - 1+o( ' 1 " , where pk can be calculated numerically. 



1 Introduction 

A triangulation is a decomposition into triangles. In this paper we are interested in triangulations 
of point sets in the Euclidean plane, where the input is a set of points in the plane, denoted V, 
and a triangulation is defined as a subdivision of the convex hull of V into triangles whose vertices 
are the points in V . 

It is a well-known fact that n points in the plane can have many different triangulations. For most 
application domains, the choice of the triangulation is important, because different triangulations 
can have different effects. For example, two important fields in which triangulations are frequently 
used are finite element methods and terrain modeling. In the first case, triangulations are used to 
subdivide a complex domain by creating a mesh of simple elements (triangles), over which a system 
of differential equations can be solved more easily. In the second case, the points in V represent 
points sampled from a terrain (thus each point has also an elevation), and the triangulation 
provides a bivariate interpolating surface, providing an elevation model of the terrain. In both 
cases, the shapes of the triangles can have serious consequences on the result. For mesh generation 
for finite element methods, the aspect ratio of the triangles is particularly important, since elements 
of large aspect ratio can lead to poorly-conditioned systems. Similarly, long and skinny triangles 
are generally not appropriate for surface interpolation because they can lead to interpolation from 
points that are too far apart. 

In most applications, the need for well-shaped triangulations is usually addressed by using the 
Delaunay triangulation. The Delaunay triangulation of a point set V is defined as a triangulation 
where the vertices are the points in V and the circumcircle of each triangle (that is, the circle 
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Fig. 1: Left: A Delaunay triangulation (k = 0). Center: an order-1 triangulation (with useful-1, 
non-Delaunay, edges in gray). Right: an order-2 triangulation, with order-1 triangles in 
light gray and order-2 triangles in medium gray. 

defined by the three vertices of each triangle) contains no other point from V . The Delaunay 
triangulation has many known properties that make it the most widely-used triangulation. In 
particular, there are several efficient and relatively simple algorithms to compute it, and its tri- 
angles are considered well-shaped. This is because it maximizes the minimum angle among all 
triangle angles, which implies that its angles are — in a sense — as large as possible. Moreover, 
when the points are in general position (that is, when no four points are cocircular and no three 
points are collinear) it is uniquely defined. However, this last property can become an important 
limitation if the Delaunay triangulation is suboptimal with respect to other criteria, independent 
of the shape of its triangles, as it is often the case in applications. 

To overcome this limitation, Gudmundsson et al. proposed higher order Delaunay triangula- 
tions [10]. They are a natural generalization of the Delaunay triangulation that provides well- 
shaped triangles, but at the same time, flexibility to optimize some extra criterion. They are 
defined by allowing up to k points inside the circumcircles of the triangles (see Figure 1). For 
k = 0, each point set in general position has only one higher order Delaunay triangulation, equal 
to the Delaunay triangulation. As the parameter k is increased, more points inside the circum- 
circles imply a reduction of the shape quality of the triangles, but also an increase in the number 
of triangulations that are considered. This last aspect makes the optimization of extra criteria 
possible, thus providing triangulations that are a compromise between well-shaped triangles and 
optimality with respect to other criteria. 

Therefore the importance of higher order Delaunay triangulations lies in multi-criteria triangula- 
tions. Their major contribution is providing a way to optimize over a — hopefully large — class of 
well-shaped triangulations. 

A particularly important subclass of higher order Delaunay triangulations are the first order 
Delaunay triangulations, that is, when k = 1. It has been observed that already for k = 1, a 
point set with n points can have an exponential number of different triangulations [16]. This, 
together with the fact that for k = 1 the shape of the triangles is as close as possible to the shape 
of the Delaunay triangles (while allowing more than one triangulation to choose from) , make first 
order Delaunay triangulations especially interesting. In fact, first order Delaunay triangulations 
have been shown to have a special structure that facilitates the optimization of many criteria [10]. 
For example, it has been shown that many criteria related to measures of single triangles, as 
well as some other relevant parameters like the number of local minima, can be optimized in 
O(nlogn) time for k — 1. In a recent paper [17], Van Kreveld et al. studied several types of more 
complex optimization problems, constrained to first order Delaunay triangulations. They showed 
that many other criteria can be also optimized efficiently for k = 1, making first order Delaunay 
triangulations even more appealing for practical use. 

For larger values of fc, fewer results are known. The special structure of first order Delaunay 
triangulations is not present anymore, which complicates exact optimization algorithms. Several 
heuristics and experimental results have been presented for optimization problems related to ter- 
rain modeling, showing that very small values of k (k = 1, . . . , 8) are enough to achieve important 
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improvements for several terrain criteria [4,5,8]. 

However, despite the importance given to finding algorithms to optimize over higher order Delau- 
nay triangulations, it has never been studied before how many higher order Dclaunay triangulations 
there can be in the hrst place. In other words, it is not known what the minimum and maximum 
number of different triangulations are, as functions of fc and n, not even for the simpler (but — in 
practice — most important) case of first order Dclaunay triangulations. 

The problem of determining bounds on the number of higher order Dclaunay triangulations is of 
both theoretical and practical interest. 

From a theoretical point of view, determining how many triangulations a point set has is one of 
the most intriguing problems in combinatorial geometry, and has received a lot of attention (e.g. 
[2, 14, 15]). Higher order Delaunay triangulations are a natural and simple generalization of the 
Dclaunay triangulation, hence the impact of such generalization on the number of triangulations 
is worth studying. 

From a more practical point of view, knowing the number of triangulations for a given fc gives an 
idea of how large the solution space is when optimizing over this class of well-shaped triangulations. 
Ideally, one expects to have many different triangulations to choose from, in order to find one that 
is good with respect to other criteria. 

Up to now, only trivial bounds were known: every point set has at least one order-fc Delaunay 
triangulation, for any k (equal to the Delaunay triangulation), and there are point sets of size n 
that have 2 e (") triangulations, already for k = 1. In this paper we present the first non-trivial 
bounds on the number of higher order Delaunay triangulations. Given the practical motivation 
mentioned above, we are mostly interested in results that have practical implications for the use 
of higher order Delaunay triangulations. Thus low values of k are our main concern. Our ultimate 
goal — achieved partially in this paper — is to determine to what extent the class of higher order 
Dclaunay triangulations (for small values of fc, which has the best triangle-shape properties), also 
provides a large number of triangulations to choose from. 

Results We study lower and upper bounds on the number of higher order Delaunay triangula- 
tions, as well as the expected number of order- A: Delaunay triangulations for uniformly distributed 
points. Let Tk(n) denote the maximum number of order-/,: Delaunay triangulations that a set with 
n points can have, and let tk(n) denote the minimum number of order-fc Delaunay triangulations 
that a set with n points can have. First we show that the lower bound tk(n) > 1 is tight. In 
other words, there are arbitrarily large point sets that have a single higher order Delaunay trian- 
gulation, even for large values of k. Next we show that, for first order Delaunay triangulations, 
Ti(n) = 2 ra_3 . Since these extreme cases do not describe an average situation when higher order 
Dclaunay triangulations are used, we then study the number of higher order Delaunay triangula- 
tions for a uniformly distributed point set. Let Rk denote the number of order-fc (and not order-z, 
for any i < k) Delaunay triangulations of a uniformly distributed point set of size n. We show 
that E[i?i] > 2' 5l "( 1+ °W), where pi is an analytically defined constant {p 1 ps 0.525785). We also 
prove that, for constant values of k, E[R k ] > 2^™( 1+ °W) , where Pk can be calculated numerically 
(asymptotics are with respect to n). The result has interesting practical consequences, since it im- 
plies that it is reasonable to expect an exponential number of higher order Dclaunay triangulations 
for any k > 1. 

Related work As mentioned earlier, there is no previous work on counting higher order Delaunay 
triangulations. A related concept, the higher order Delaunay graph, has been studied by Abellanas 
et al. [1]. The order-k Delaunay graph of a set of points V is a graph with vertex set V and an 
edge between two points p, q when there is a circle through p and q containing at most k other 
points from V . Abellanas et al. presented upper and lower bounds on the number of edges of this 
graph. However, since a triangulation that is a subset of the order-fc Delaunay graph does not 
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need to be an order-/. - Dclaunay triangulation, it is difficult to derive good bounds for higher order 
Delaunay triangulations based on them. 

There is an ample body of literature on the more general problem of counting all triangulations. 
Lower and upper bounds on the number of triangulations that n points can have have been 
improved many times over the years, with the current best ones establishing that there are point 
sets that have 0(8.48") [2] triangulations, whereas no point set can have more than 0(43") [15]. 

In relation to our expected case analysis of the number of higher order Delaunay triangulations, it 
is worth mentioning that many properties of the Delaunay triangulation — and related proximity 
graphs — of random points have been studied. The expected behavior of properties of the Dclaunay 
triangulation that have been considered include the average and maximum edge length [3, 12], the 
minimum and maximum angles [3], and its expected weight [6]. Expected properties of other 
proximity graphs, such as the Gabriel graph and some relatives, are investigated in [7,9, 11]. 

Outline This paper is structured as follows. The next section presents some previous results 
related to higher order Delaunay triangulations, needed for the following sections. In Section 3 we 
give lower and upper bounds for the number of higher order Delaunay triangulations. Section 4 
deals with the expected number of higher order Delaunay triangulations. Finally, some concluding 
remarks are made in Section 5. 

2 Higher order Delaunay triangulations 

We begin by introducing higher order Delaunay triangulations more formally, and presenting a 
few properties that will be used throughout the paper. From now on, we assume that point sets 
are in general position. 

Definition 1: (from [10]) A triangle Auvw in a point set V is order-k Delaunay if its circumcircle 
C(u,v,w) contains at most fc points of V . A triangulation of V is order-k Delaunay if every 
triangle of the triangulation is order- k. 

Note that if a triangle or triangulation is order- A:, it is also order-/;' for any k! > k. A simple 
corollary of this is that, for any point set and any k > 0, the Delaunay triangulation is an order-fc 
Delaunay triangulation. 

Definition 2: (from [10]) An edge uv is an order-fc Delaunay edge if there exists a circle through u 
and v that has at most k points of V inside. An edge uv is a useful order-fc Dclaunay edge (or 
simply useful-k edge) if there is an order-/.' Dclaunay triangulation that contains uv. 

The useful order of an edge can be checked using the following lemma, illustrated in Figure 2. 

Lemma 1: (from [10]) Let uv be an order-fc Delaunay edge, let s\ be the point to the left 1 of v&, 
such that the circle C{u,v,s\) contains no points to the left of vu. Let S2 be defined similarly 
but to the right of vu. Edge uv is useful-fc if and only if Auvsi and Auv.S2 are order-fc Delaunay 
triangles. 

The concept of a fixed edge is important in order to study the structure of higher order Dclaunay 
triangulations. 

1 We sometimes treat edges as directed, to be able to refer to the right or left side of the edge. The left side of 
vu denotes the halfplane defined by the line supporting IS, such that a polygonal line defined by v, u and a point 
interior to that halfplane, makes a counterclockwise turn. In the right side, the turn is clockwise. 
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Fig. 2: The useful order of edge uv is determined by the lowest order of triangles Auvsi and 
Auvs-2- In the example, the (lowest) useful order of uv is max{3, 1} = 3. 

Definition 3: Let V be a point set and T its Dclaunay triangulation. An edge of T is k-fixed if it 
is present in every order- k Delaunay triangulation of V . 

Some simple observations derived from this are that the convex hull edges are always fc-fixcd, for 
any k, and that all the Dclaunay edges are 0-fixcd. 

First order Delaunay triangulations have a special structure. If we take all edges that are 1-fixed, 
then the resulting subdivision has only triangles and convex quadrilaterals (and an unbounded 
face). In the convex quadrilaterals, both diagonals are possible to obtain a first order Delaunay 
triangulation (see Figure 1, center). We say that both diagonals are flippable, and similarly we 
call the quadrilateral flippable. More formally, based on results in [10], we can make the following 
observation. 

Observation 1: Let e be a useful ordcr-1 Dclaunay edge in an order-1 Delaunay triangulation, such 
that e is not a Delaunay edge. Then flipping e results in a Delaunay edge. Moreover, the four 
edges (different from e) that bound the two triangles adjacent to e are 1-fixed edges. 

An implication of this special structure is that instead of counting triangulations, we can count 
flippable quadrilaterals or, cquivalently, useful-1 edges that are not Dclaunay. 

Corollary 1: Let V be a point set. If V has q flippable quadrilaterals, then V has exactly 2 q order-1 
Delaunay triangulations. 

For k > 1, the structure is not so simple anymore and it seems difficult to provide an exact 
expression for the number of order-fc Delaunay triangulations in terms of the number of useful-fc 
edges. However, we can derive a lower bound by combining a number of known results, as follows. 
First we need some extra definitions and previous results: 

Definition 4: (from [10]) The hull of an order-fc Delaunay edge uv (k > 1) is the closure of the 
union of all Delaunay triangles whose interior intersects uv. 

Lemma 2: (from [10]) The hull of an order-fc Dclaunay edge (k > 1) is a simple polygon consisting 
of at most 2k + 2 vertices. 

Lemma 3: (from [10]) Let uv be a useful-/,- edge (and not useful-? for any i < k), with k > 1. There 
exists an order-fc (and not order-i for any i < k) Delaunay triangulation of the hull of uv that 
contains uv. 

Lemma 4: Let uv be an order-0 edge. The number of useful-fc edges (fc > 1) that intersect uv is 
at most (2k + l) 2 . 
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Proof: It follows directly from the proof of Lemma 8 in [10]. Kl 

We have now the necessary tools to prove the following lower bound on the number of order-fc 
triangulations, expressed as a function of the number of useful-/e edges. 

Lemma 5: Let V be a point set and let ek, for k > 1, be the number of useful-/,- edges (which arc 
not useful-i for any i < fc) of V . Then V has at least 2 Ck ' Ck — 1 order-fc (and not order-z, for any 
i < k) Delaunay triangulations, where Ck — (4fc + l)(2fc + l) 2 + 1. 

Proof: Let Ek denote the set of useful-fc edges (which are not useful-i for any i < k) of V and let 
efc denote the cardinal number of Ek . We select a subset E' k of the edges of Ek in the following 
way: We pick an edge e of Ek. we remove all the edges in Ek whose hull intersects the hull of e in 
at least one Delaunay triangle, and we repeat until Ek does not contain any edge. 

Let e' be an edge in E' k and e be an edge in Ek whose hull intersects the hull of e' in at least one 
Delaunay triangle T. Then e intersects at least one edge of T. By Lemma 2, all Delaunay triangles 
included in the hull of e' contain at most 4fc + 1 edges. By Lemma 4, each of them intersects at 
most (2fc + l) 2 useful-fc edges. Hence if e' is selected, at most (4k + l)(2k + l) 2 edges in Ek are 
removed. Therefore E' k contains at least efc/((4fc + l)(2fc + l) 2 + 1) edges. 

Each non-empty subset of E' k gives rise to a different order-fc (and not order-i, for any i < fc) 
Delaunay triangulation proceeding as follows: If an edge e is in the subset, we triangulate the 
hull of e as in Lemma 3, that is, we use an order-fc (and not order-z, for any i < fc) Delaunay 
triangulation containing e. If an edge e is not in the subset E' k , we triangulate the hull of e using 
the Delaunay triangles crossed by e. Finally, we complete the triangulation by adding Delaunay 
triangles in the regions that have not been triangulated (that is, computing a constrained Delaunay 
triangulation). This construction is consistent because the hulls of the edges in E' k can only 
intersect in points and boundary edges, and because the boundary edges of the hulls belong to 
the Delaunay triangulation. Kl 



3 Lower and upper bounds 

In this section we derive upper and lower bounds on the number of first order Delaunay triangu- 
lations. As mentioned in the introduction, due to the practical motivation of this work, we are 
mostly interested in lower bounds. However, for completeness and because the theoretical question 
is also interesting, this section also includes a result on upper bounds. 

The main question that we address in this section is: what is the minimum number of higher 
order Delaunay triangulations that /?. points can have? Are there arbitrarily large point sets that 
have only 0(1) higher order Delaunay triangulations? To our surprise, the answer to the second 
question is affirmative. 

The lemma below presents a construction that has only one higher order Delaunay triangulation, 
regardless of the value of fc, for any fc < [n/3j — 1. Note that this implies that for any value of fc 
of practical interest, there are point sets that have no other order-fc Delaunay triangulation than 
the Delaunay triangulation. 

Lemma 6: Given any n > 6 and any fc such that fc < [n/3j — 1, there are point sets with n points 
in general position that have only one order-fc Delaunay triangulation. 

Proof: We give a construction with n points that can be shown to have only one order-fc Delaunay 
triangulation, for any value of fc < [n/3j — 1. The construction is illustrated in Figure 3. 

To simplify the explanation, in the following we assume that ?? is a multiple of 3. Since any order-?' 
Delaunay triangulation is also order-fc for all fc > ?', for the proof it is enough to use k = n/3 — 1. 
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Fig. 3: Construction of a point set (left) whose only order-fc Delaunay triangulation is the Delaunay 
triangulation (right). 

We start with a triangle AS1S2S3. Then we add three groups of points, which we will denote with 
letters p, q, and r. The points in the first group are denoted p\ . . .p m , where m = n/3. These 
points are initially placed on a circle C p that goes through S3, as shown in the figure; they are 
sorted from top to bottom. The second group comprises k points q\ . . .qk, placed very close to 
each other on a circle C q through S3, as shown in the figure. In addition, we must also make 
sure that the points q\ . . . qu are close enough to S3 in order to be contained inside C(si, S2,Pi)- 
Finally, the points in the third group, r± . . . rk—i, are placed very close to each other on a third 
circle C r , which goes through 82- The important properties of these circles are: (i) C p contains 
S2 and all the points rj, (ii) C r and C q contain all the points of the type pt, and (iii) C(si,S2,Pi) 
contains S3 and all points of type qi. 

Clearly, the point set as constructed is degenerate, but this can be easily solved by applying a 
slight perturbation to each point, without affecting the properties just mentioned. Moreover, the 
perturbation can be made such that the Delaunay triangulation of the point set looks like the one 
in the right of Figure 3. 

We now argue that all the edges in the Delaunay triangulation are fc-fixed, by considering the 
different types of edges that, potentially, could cross a Delaunay edge to make it non-fixed. Suppose 
an edge of the shape R\p\ is not fc-fixed. Then there must be some triangulation in which the edge 
is crossed by some other useful order-fc edge. Such edge can be of three types: (i) it connects two 
points pj, pk, (ii) it connects two points pj, qk (or S3), or (iii) it connects two points pj, rj. (or 
.52)- An edge of the type pjpk that crosses sTpJ must be an edge of the shape PiPi+2 or force a 
such an edge to appear in the triangulation. However, the circumcircle of the triangle defined by 
any three consecutive points pi, p%+i, Pi+2 contains at least k + 1 points because it is a slightly 
perturbed version of C p . Thus no such edge can be part of an order-fc triangulation. A similar 
situation occurs with any edge of the shape pjqk, since it forces a triangle of the form Apipiq m 
(or Apipis^). Finally, edges of type pfnl force a triangle of the form Asip^r; (or AsiPiS2), whose 
circumcircle includes at least as many points as contained in C(s\, S2,Pi), hence cannot be part 
of an order-fc triangulation either. Therefore all the edges of the shape sTJjJ are fc-fixed. Similar 
arguments can be used to show that the edges in the other groups are also fc-fixed, hence no other 
order-fc triangulation can exist. S3 

Having determined that some point sets can have as little as one first order Delaunay triangulation, 
it is reasonable to ask what is the maximum number of first order Delaunay triangulations that 
a point set can have. The following lemma gives a precise — and tight — bound on the maximum 
number of first order Delaunay triangulations. 
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Fig. 4: Construction achieving the maximum number of first order Delaunay triangulations. Left 
and right: point set and fiippable quadrilaterals, for points not in general position. 

Lemma 7: Every point set V with n points in general position has at most 2™~ 3 first order Delaunay 
triangulations, and this bound is tight. 



Proof: To see that no point set can have more than n — 3 fiippable quadrilaterals, observe that 
the subdivision of the convex hull of V induced by the fixed edges is a plane graph. It follows from 
Euler's equation that any triangulation has at most 2n — 5 triangles. Since each quadrilateral is 
formed by two triangles, there can be at most n — 3 quadrilaterals. 

Now we give a construction with n (for n any multiple of 4) points that has n — 3 fiippable quadri- 
laterals, thus a total of 2"~ 3 first order Delaunay triangulations. The construction is illustrated 
in Figure 4, and consists of a series of points placed on the vertices of concentric squares with 
the same orientation. Clearly, the edges in Figure 4 are Delaunay edges and the four vertices 
of each quadrilateral are cocircular. If we apply a small perturbation to the point set so that it 
reaches a general position, one of the diagonals of each quadrilateral becomes a Delaunay edge, 
while the other one becomes a useful order-1 (non-Delaunay) edge. Therefore, all quadrilaterals 
are fiippable. 



4 Expected number of triangulations 

Let V be a set of n points uniformly distributed in the unit square. In this section we give lower 
bounds on the expected number of higher order Delaunay triangulations of V . 

Note that the events that four points in V are cocircular and that three points in V form a right 
angle happen with probability zero, and hence we can safely ignore these cases. Throughout this 
section we will use the notation x ~ y if x = y(l + o(l)). 

We start with first order Delaunay triangulations. We aim to compute the probability that two 
randomly chosen points u,v in V form a useful-1, non-Delaunay, edge. Assume that the edge is 
directed vi,. Let w be the point to the left of vu, such that the circle C(u, v, w) contains no points 
to the left of vu, and let t be the point to the right of vu, such that the circle C(u,v,t) contains 
no points to the right of vu. Let £ be the event defined as follows: the edge uv is useful-1 (but 
not Delaunay), t is to the right of vu and the circle C(u, v, t) contains no points of V to the right 
of vu, w is to the left of vu and the circle C(u, v, w) contains no points of V to the left of vu". It is 
well-known that uv belongs to the Delaunay triangulation of V if and only if Zuwv + Zutv < tt. 
Thus the event £ can be decomposed into the disjoint union £ = £\ U<?2 U^3, where £\ denotes the 
event £ with the additional conditions that Z.uwv > ir/2, Autv > tt/2, £2 denotes the event £ with 
the conditions that Zuwv < n/2, Zutv > n/2, and £3 denotes the event £ with the conditions 
that Zuwv > 7r/2, Zutv < n/2 (in all cases we must have Zuwv + Zutv > tt). Consequently, 
V[£]=P[£ 1 }+P[£ 2 }+V[£ 3 }. 



Lemma 8: P[£i] - a/n 3 and P[£ 2 ] - c 2 /n 3 , where a = 0.23807 and c 2 = 0.40675. 



Fig. 5: Left: the event £\ and the regions A w , B w , A t and B t . Middle: in the event £\, the region 
A w is a circular sector minus a triangle. Right: the event £2 and the regions A w , B w , Af 
and B t . 



Proof: Let us first compute P[£i]. 

Let A w (respectively, A t ) be the interior of the set consisting of all points in C(u, v, w) DC(u, v, t) 
that are to the left (resp. right) of vu. Let B w (respectively, B t ) denote the interior of the set 
containing all points in C(u,v,w) (resp. C(u,v,t)) that are to the right (resp. left) of vu and do 
not lie in A t (resp. A w ) (see Figure 5, left). Since w is the point such that the circle C(u,v,w) 
contains no points to the left of vu, the region A w is empty of points in V . In order for the edge 
uv to be useful- 1, the region B t also has to be empty of points in V . Analogously, under the 
hypothesis of £\, the regions A t and B w are empty of points in V . It is not difficult to see that 
the reverse implications also hold. Therefore, the event £\ is equivalent to the event that A w , Bt-, 
At, and B w do not contain any point in V . 

Now let us denote by r w the radius of the circle C(u, v, w) and by £ the length of the edge uv (see 
Figure 5, center). A straightforward calculation leads to the following expression for the area of 
A w : 



area(j4„,) = r w arcsin 




In order to compute ¥[£], we will be interested in having certain areas being empty of n points, 
which happens with probability (1 — area {A)) n (for A the area in question). Since the contribution 
of areas A of size 0(1) is 0(X n ) for some < A < 1 (which is far less than the asymptotic value 
of the integrals, as we shall see below), for any constant j we can safely assume in the integrals 
below the asymptotic equivalence (1 — area (.A)) n_J ' ~ e — arca(A)^ w jthout affecting the first order 
terms of the asymptotic behavior of the integral 2 . 

Observe that t may take values from to y/2 and that the probability density of the event 
|m;| = £ is 2iztdt. Notice also that r w € (£/2,+oo) and that the event of having a radius r w 

has probability density I —2r w arcsin ( J H — , o ) dr w , since it corresponds to the negative 



derivative — f (r w )dr w of the function f(r) = r 2 arcsin (^) — | y r 2 — ^. 

Denoting by r f the radius of the circle C(u,v,t), we obtain analogous expressions for r t . 

Now we have all the necessary ingredients to develop an expression for P[£i] : 



2 In fact, this formula arises in a homogeneous Poisson point process of intensity n in the unit square, and it 
is not surprising that both distributions give the same asymptotic results (see the ideas of depoissonization given 
in [13]). 
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2^ 



•[Sx] ~ / 2tt£ / / -2r w arcsin 

JO J 1/2 J 1/2 

t\ . £r t \ (1) 
-2r t arcsin | - — 



2? 



e -n(irrl +7rr t 2 -r I arcsin( ^ ) + § \j - ^ -r f 2 arcsin( ^ )+ f \jr\-^) ^ , ^ 

Since classical methods for asymptotic integration fail for the integral given by (1) (the derivative 
of the exponent is infinity at the point where the exponent maximizes), we apply the following 
change of variables: £/2 — a/ 'y/n, r t — bj \/n, r w — cj •Jn. The integral (1) then becomes (replacing 
the integration limit \[2~n by oo, which can be done since the dominant contribution comes from 
small values of a) 




Given that it does not seem possible to evaluate this integral symbolically, we resort to applying 
numerical methods. For reasons of numerical stability (especially in the case of ¥[£2] below) we 
apply another change of variables b — sin (°y 2 ) > c = sin(6> /2) ' 



mr „ n 1 f°° f w r n fa 2 9cot(0/2)~2a 2 \ f a 2 <j cotia 12) ~ 2a? 



n 3 Jo Jo Jo V 1 — cos(0) J \ 1 — cos(cr) 

^ 2 (T^-cot(e/2) + T ^- T -cot( ff /2)) d(j ^ a ^ 



Solving this integral numerically (using Mathcmatica) , we obtain that P[£i] ~ ci/n 3 , where c\ — 
0.23807. 

Let us now consider £ 2 . Let A w , B w , A t and B t be defined as in the event £\ (see Figure 5, right). 
By the same arguments, the event £ 2 is equivalent to the event that the regions A w , B tl A t , and 
B w are empty of points in V . Analogous observations as in the previous case yield 



ry/2 poo poo 

'[£2] ~ / 2tt£ / / I 2r w ix — 2r w arcsin 

JO J 1/2 Jr m 



t \ £r t 
-2r t arcsin I — ^ 
. 2r t / 



t \ £r„ 



' -hi 4 



-n(7rr 2 -r 2 arcsin ( 277) + ! y^? ~ T + r ™ arcsin( ji-) - 1 yVj, - ^ ) 



dr t dr w d£. (2) 



As before, we apply the substitution £/2 = a/^/n, r t — b/y/n, r w = cj 1 \fn to the integral (2) and 
obtain 
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- 7r b 2 +b 2 arcsin ( f ) - a Vb 2 -a 2 - c 2 arcsin ( § ) +a Vc 2 -a 2 ^ ^ c dfl 



For reasons of numerical stability, we again apply the change of variables b — sin ^°y 2 ) an< ^ c 
and obtain 



i(9/2) 



1 f°° f n f e a fa 2 {9 - 2tt) cot(0/2) - 2a 2 \ /a 2 a cot(a/2) - 2a g 



o Jo Jo 



1 — cos(0) / \ 1 — cos(tr 



a-(cot(e/2)- T= ^| iW + r 2^ T - co t(<7/2)) da . d0da _ 



Solving this integral numerically (using Mathcmatica) , we obtain that ¥[£2] ~ c'2/n 3 , where c 2 = 
0.40675. IE 

Denote by £7i the random variable counting the number of useful- 1 (and not Delaunay) edges. We 
have the following corollary: 



Corollary 2: E[Ui] ~ c -^^ =: dm, where di = 0.525785. 



Proof: Since £2 and £3 are symmetric, we obviously have that P[£?,} = ¥[£2] - Hence, P[£] ~ Cl ^ 3 C2 . 
Since for a fixed edge uv there are (n — 2) (n — 3) ~ n 2 ways to choose the points w and t to the left 
and to the right of vu, and these events are all disjoint, the edge uv is useful- 1 (and not Delaunay) 
with probability £l± ^. Hence, E[C/i] - Q^^r 2 - ~ £3 ^r a = : ^n, where d x = 0.525785. M 

Recall that Rk denotes the number of order-fc (and not order-i, for any i < k) Delaunay triangu- 
lations of a uniformly distributed point set. We can now state the following theorem: 

Theorem 1: Given n points distributed uniformly at random in the unit square, E[i?i] > 2 Pl "( 1+ °W, 
where pi = 0.525785. 



Proof: By Corollary 1, E[i?i] = E[2 Ul \. Now, by Jensen's inequality, E[2 Ul ] > 2^ Ul \ and the 
result follows by Corollary 2. H 

Combining the ideas for the case k — 1 with the result from Lemma 5, we obtain the following 
generalization for constant values of k: 



Theorem 2: Given n points distributed uniformly at random in the unit square, for any constant 
value of k, E[R k ] > 2^' l ( 1 +°( 1 », where Pk is a constant that can be calculated numerically. 



Proof: Denote by U k the number of useful-fc edges (which are not useful-i for any i < k) for any 
constant k > 1. We want to know the value of E [[/*;]. 

In order for an edge uv to be useful-fc (and not useful- i for i = 0, . . . , k — 1), using the notation of 
Figure 5, first observe that the regions A w and A t have to be empty of points. Moreover, either 
the region B t \ A w has to contain exactly k — 1 points (w is excluded) , whereas the region B w \ A t 
can contain any number of points i = 0, . . . , k — 1 (t is excluded), or vice versa. For any constant 
i, the probability of having exactly i points in an area A of size o(l) (as before, for constant i only 
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such areas count for the asymptotic behaviour of the integrals) is ~ e (nA) ! /i!. Thus, defining 
the events £1, £2, and £3 analogously as in Section 4, 



^ 1 f-^Pl poo />oo 

\Ei] ~ ^2 / 2?r W / I - 2r ™ arcsin 



-Zr t arcsm 



2r„ 



2r 



e -n(.irr% +^-rl arcsin( ^- ) + § ^jr^-^-rl arcsin( ^ ) + f \jr\-^) 

(area (5 t \ A») ™) fe_1 (area \ A t ) n) 1 - — l ——dr t dr w d£ 

(k — l)h\ 



+ ^ / 2tt£ / / I — 2r w arcsin 



^ J Jrl, - f 



-2r t arcsin 



2r t J l~2 1 



^-nK+irJ-rJ, arcsin( ^- ) + § - -r t 2 arcsin ( ^ ) + § \Jr^-^-) 



(area (£ t \ A w ) n)'(area (B w \ A t ) n) k 1 - - TT -dr t dr w d£. 

(k — l)h\ 



Now, since 



i \ £ / 2 i 2 2 / e \ £ / „ p 

area (5 4 \ A„) = r t tt - r t arcsin ( — I + - y r; ?- — - r w arcsin I — I + - ^ r 2 w - — 



and 



£ / _ £ 2 2 / M * / 2 ^ 2 

area(B lu \ A t ) = r^7r - r.^ arcsin - — + -\/r^ — r t arcsin — + -W rf 



2r w J 2 V tu 4 c V2r f y 2V 1 4 

after applying the substitutions £/2 = a/^/n, rt = b/^/n, r w = c/^/ri 1 in these new factors n 
disappears and the integral again yields 0(l/n 3 ). The same argument also holds for P[£2], and 
by the same reasoning as in the case of useful-1 edges, we obtain that the expected number of 
uscful-A: edges (that are not useful-i for any i < k) is d k n for any constant fc. We point out that 
using our formula the constant d k can be calculated numerically. 

By Lemma 5, R k > 2 U ^ C " - 1, where C k = {ik + l)(2k+ 1) 2 + 1. Therefore E[R k ] > E[2 u ^ Ck ] - 1 
and as before, by Jensen's inequality, E[2 u "/ Ck ] > 2 E \- U »V C » . H 

5 Discussion and further work 

We have given the first non-trivial bounds on the number of higher order Delaunay triangulations. 
We showed that there are sets of n points that have only one higher order Delaunay triangulation 
for values of k < [n/3\ — 1, and that no point set can have more than 2"~ 3 first order Delaunay 
triangulations. Moreover, we showed that for any constant value of k (in particular, already for 
k — 1) the expected number of order-/,: triangulations of n points distributed uniformly at random 
is exponential. This supports the use of higher order Delaunay triangulations for small values of k, 
which had already been shown to be useful in several applications related to terrain modeling [8] . 
From a more theoretical perspective, it would be interesting to obtain tighter bounds on the 
expected number of order-A; Delaunay triangulations for uniformly distributed points. 
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